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ABSTRACT 

We analyzed vibrational stability of metal-poor low-mass main-sequence stars due to 
the e-mechanism. Since outer convection zones of the metal-poor stars is limited only 
to the very outer layers, the uncertainty in treatment of convection does not affect 
the result significantly. We found that the dipole gi- and g2-modes certainly become 
unstable due to the e-mechanism for Z < 6 x 10^'*. Besides that, we found that as the 
metallicity decreases the mass range of the e-niechanism instability extends toward 
higher mass. 
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1 INTRODUCTION 
T— I ■ 

^ ■ Structure and evolution of stars born in metal-poor environ- 
] ments or when heavy elements were siginificantly deficient 
, in the Universe are considerably different than for stars born 

■ much later with higher heavy-element abundances. This is 
[ mainly because opacities and nuclear reactions are highly 

■ dependent on metallicity. Opacity decreases with metallic- 
f"*) ■ ity because of lack of absorption by heavy elements. Hence 
fS| . it makes luminosity of the star higher. To maintain energy 
T-H ■ equilibrium in this situation, the stars with low metallicity 

' need to be compact compared with the metal-rich stars. As 
. ^ , a result, the main-sequence of metal-poor stars on the HR 

■ diagram moves toward the bluer and higher luminosity side. 
$H ] With this, the convective envelope of stars with ~ IMq be- 
5^ ■ comes t hiner a nd limited to the very outer layer close to the 

surface (ISonoi fc Shibahashi 2011j, hereafter Paper I). 

The decrease in metallicity also makes CNO-cycle en- 
ergy generation less efficient. In the Population I case, the 
pp-chain is dominant for < 1.2Mq, while the CNO-cycle 
is infiuential for the more massive stars. But as the metal- 
licity and the CNO abundance decreases, the CNO-cycle 
becomes replaced by the pp-chain for more massive stars. 
In the Population III case, the nuclear energy source is only 
the pp-chain at the ZAM S stage for stars with < 13Mq 
l|Sonoi fc Shibahashill2012l . hereafter Paper II). 

In the case of the pp-chain burning, we should note 
the vibrational instability of low-degree low-order g modes 
due to the e-mechanism. Such instability and the resul- 
tant material mixing was once prop osed as a possible so - 
lution to the solar neutrino problem (jPilke fc Goughlll972h . 
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and the detailed numerical calculations of linear stability 
analyses demonstrated that such instability is likely to oc- 
cur in a certain ea rly evolutionary stage of the Sun and 
the solar- like stars (Christcnsen-Dalsgaard, Dilko & Gou^ 
11974; Bo urv et al.l ll975: Shibahashi ct al. 1975 : Nool s et akl 
[l976l ). The presence of a convective envelope, which occupies 
the outer 20-30% of the stellar radius, however, has made 
it hard to reach a definite conclusion on the vibrational sta- 
bility because of uncertainty in treatment of the convective 
envelope. 

However, the situation is different in metal-poor stars. 
As described above, with decreased metallicity, convective 
envelopes of the solar-like stars become thin enough that the 
uncertainty in treatment of convection does not significantly 
affect results of the vibrational stability analyses. 

Indeed, the present authors showed that metal-free Pop- 
ulation HI stars are vibrationally unstable against dipole 
low-order g modes due to the e-mechanism of the pp-chain 
(Papers I and II) . To extend our analyses to stars with Z 
but still having only very thin convective envelope, we ex- 
amined the stability of stars with low metallicity by a non- 
adiabatic analysis, and determined the upper limit metal- 
licity for which this instability appears without the uncer- 
tainty. 



2 EVOLUTIONARY EQUILIBRIUM MODELS 

We adopted the same code as in Paper II, MESA 
l|Paxton et al.ll201lf ). to calculate stellar evolution. Evolu- 
tionary models were constructed by calculating from the 
pre-main-sequence Hayashi phase with the initial hydro- 
gen abundance, Xq = 0.75, and the metallicities, Z — 
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Figure 1. Top and bottom boundaries of an outer convection 
zone for ZAMS models with Z = 4x 10"*, 1 x IQ-^ and 2x 10^^ 
The initial hydrogen abundance is set to be Xq = 0.75. (^4 color 
version of this figure is available in the online journal). 



Figure 3. The fraction of nuclear energy by the CNO-cycle to 
the whole nuclear energy at the stellar center for 1.O-2.4M0 with 
a O.2M0 step. The abscissa is the central hydrogen abundance 
Xc, which is an indicator of stellar evolution. This range is set to 
be from X^ = 0.75 to 0.20. 
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Figure 2. Evolutionary tracks for 0.9-1.1Mq with Xo = 0.75, 
Z = 1, 2, 4, 6, 8 X lO"** and 1 X 10~^. For evolutionary models on 
the thick parts of the tracks, the convective envelope is limited to 
the outer layer above r/R = 0.99. The thick dashed line means 
the red edge of the region where such models are located (^4 color 
version of this figure is available in the online journal). 



1,2,4,6,8x10"" and 1x10"^ The mixing length parameter 
was set to be omlt = 1-0. 

Fig.[T]shows top and bottom boundaries of an outer con- 
vection zone for ZAMS models with different metallicity. As 
described in Sec. [T] with decreasing metallicity, outer con- 
vection zones for ~ IMq stars become thinner, limited to 
outer layers above r/R ~ 0.99, and as thin as or thinner than 
of population I stars with 1.6-2.0Mq, which correspond to 
the blue edge of the S Scuti instability strip. The vibra- 
tional stability of those stars has been well-dete rmined with- 
out uncertainty in treatment of convectio n fe.g. IStellingwer j 
1 19791 : iTsvetkovlfigj^ : |Pamvatnvkhll2000l ). In this study, we 
restricted our stability analysis to evolutionary models hav- 
ing very thin convective envelopes which are confined to the 
outer layers above r/R — 0.99. Fig. [2] shows evolutionary 
tracks for 0.9-1. IMq, and the thick parts of the tracks cor- 



respond to such models. The thick dashed line means the 
red edge of the region where such models are located. 

On the other hand, properties of nuclear reactions also 
change with metallicity. Fig. [3] compares contribution of the 
CNO-cycle to the whole nuclear energy generation at the 
stellar center between Z = 1 x 10~^ and 1 x 10""*. The 
contribution of the CNO-cycle decreases with the metallicity 
for the same mass star. Energy generated by the pp-chain 
occupies more than 80% for < 2.2M0 with Z = 1 x 10"* at 
the ZAMS stage, while for < 1.6Mq with Z = lx 10"^. 



3 LINEAR FULLY NONADIABATIC 
STABILITY ANALYSIS 

We performed a linear fully nonadiabatic analysis for the 
above equilibrium models by usi ng a Henyey-type code de- 
veloped in Paper II. By following lUnno et all l| 19891 ). we lin- 
earized the equations of continuity, motion, energy conser- 
vation and radiative diffusion, and the Poisson's equation, 
while any perturbation terms are expressed in terms of a 
combination of a spatial function and a time varying func- 
tion. The latter is expressed by exp(io't), where a denotes 
the eigenfrequency. The former, the spatial part, is decom- 
posed into a spherical harmonic function, which is a func- 
tion of the colatitude and the azimuthal angle, and a radial 
function. Equations governing the radial functions lead to a 
set of sixth order differential equation, of which coefficients 
are complex, including terms of frequency a. Together with 
proper boundary conditions, this set of the equations forms 
a complex eigenvalue problem with an eigenvalue a. The real 
part of a, ctr, represents the oscillation frequency, and the 
imaginary part of ct, cti, gives the growth rate or the damping 
rate, depending on its sign. We adopted the "frozen-in con- 
vection" approximation, that is, simply neglected convective 
flux perturbation. This approximation is acceptable to this 
study since outer convection zones of the chosen models are 
very thin and the result of the stability analysis may not 
depend significantly on treatment of convection. 



© 2002 RAS, MNRAS 000, [T]-?? 



e-mechanism instability of metal-poor low-mass stars 3 



o 

o 

o 



1 

0.5 

-0.5 



-0.5 



"Si 










2.0 

1 1 1 1 1 1' 1 


2.3,/'^.4 

1 1 1 


1 1 1 1 1 1 1 I'l f 




Z=lxlO''' 








1.0Mj„„ 


1 1 1 1 1, 




1.1 












1.2 

1.4 


1.8 





Table 1. Properties of the gi-mode shown in figure[5] 
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Figure 4. Variation of growth rates of the dipole{i = 1) gi- (top) 
and g2-modes (bottom) for different mass stars with Z = 1 X 10"*. 
The abscissa is same as figure [J] 



For a stability analysis relevant to the e-mechanism, 
we should adopt temperature and density dependences 
of nuclear reactions through the perturabation, et = 
{dlne/d\nT)p and Sp = {d\ne/d\n p)t, which are dif- 
ferent from those for evolutionary time-scale. For exam- 
ple, as for the pp-I branch in the pp-chain, the tem- 
pearture dependence of energy generation all through the 
branch in equilibrium is governed by the slowest reaction, 
^H(^H,e+z/e)^H, and et 4 at logT = 7. But we should 
adopt the effective temperature dependence of nuclear re- 
action through the perturbation, whi ch is mainly governed 
by ^He(^He,2^H)'^He , and £t ^ H llDilke fc Goughl ll972l: 



19891 ). We 



iBourv fc Noe3 1 19731 : lUnnd Il975l : lUnno et all 

evaluate such effective dependences of the pp-chain and the 
CN-cycle through the perturbation, separately, and then av- 
erage them with their contribution to the total nuclear en- 
ergy generation to get the net values. More details are de- 
scribed in Paper II. 



4 RESULTS 

4.1 Variation of stability with stellar evolution 
and stellar mass dependence 

In this section, we discuss variation of stability with stel- 
lar evolution and stellar mass dependence of the stability 
while showing results in the Z = 1 x 10~* case. Fig.Ushows 
variation of growth rates of the dipole {I — 1) gi- and g2- 
modes for stars with Z — I x 10~*. The positive value of 
—<Ji/{2n) means instability. Except for the lower mass stars, 
the g modes are stable at the ZAMS stage, become unsta- 
ble in the middle stage of the core hydrogen burning, and 
eventually become stable again. As well as the Population 
III case introduced in Paper I and II, the gi- and g2-modes 
become most unstable or least stable at Xc = 0.5-0.6 and 
at Xc = 0.6-0.7, respectively. This delicate change of stabil- 
ity is caused by the variation in amplitude distribution with 
stellar evolution. 

In the case of the pp-chain burning, at the stellar centre 
^He is consumed by 3He(^He,2^H)*He or by ^He(*He,7)'^Be 
immediately after it is generated. However, since the ^He 
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reactions are highly sensitive to temperature, it does not 
efficiently occur in the outer part of the nuclear burning 
core. As a consequence, ^He accumulates in an off-centered 
shell. Hence the most favorable situation for the vibrational 
instability is that the temperature perturbation has a large 
amplitude in such an off-centered '^He shell. 

Table [T] and Fig. [5] show properties of the gi-mode for 
the 1.2Mq star with Z = 1 x 10~* at different evolutionary 
stages. The core of the metal-poor low-mass stars is convec- 
tive at the ZAMS stage even with the pp-chain burning be- 
cause of the high central temperature, and the gravity waves 
are evanescent there. Note that the convective core might 
hardly affect the vibrational stability since convective time- 
scale there (~ yr) is much longer than the oscillation period 
(~ hr). In such a situation, it is pl ausible that convective 
flux does not react t o pulsation (e.g. Goldreich fc NicholsonI 
1 19771 : iPesnelll 1 19871 : iGuzik et al.l |2000| ). As the convective 
core shrinks with stellar evolution, the gravity waves be- 
come to propagate in the deep interior. The g mode then 
becomes to have a large amplitude around the '^He shell, 
and unstable due to the e-mechanism of the '^He-'^He reac- 
tion. For the higher mass stars, on the other hand, the pp-II 
and HI branches and the CNO-cycle cotribute to the nu- 
clear energy generation due to high temperature. Collisional 
reactions belonging to them have high temperature depen- 
dence comparable with that of the ^He-'^He reaction, and 
contribute to the g mode instability. In the later stage, how- 
ever, the g mode becomes to have a p mode-like behaviour in 
the envelope. Then, flux dissipation there exceeds the exci- 
tation by the e-mechanism, and the g mode becomes stable 
again. 

We found that the above instability appears for stars, 
for example, with < 2.3Af0 in the Z = 1 x 10"* case. The 
more massive stars, on the other hand, keep substantial size 
of the convective core because of the dominant contribution 
of the CNO-cycle rather than the pp-chain. In this situation, 
the gravity waves cannot propagate in the nuclear burning 
core enough for the vibrational instability. Instead, just out- 
side of the nuclear burning core, the amplitude of the g mode 
is relatively large, and strong damping is induced as shown 
in figure |6l The g modes are then not destabilized indeed in 
those stars. 

For the less massive stars, for example, ones with ^ 
I.IMq in the gi-mode case, the growth rate strongly in- 
creases in the late stage of the core hydrogen burning (figure 
|4| . Since such stars are located inside the classical instability 
strip, the g mode is excited by the K-mechanism of helium 
ionization. As the g mode becomes to have a p mode- like be- 
haviour in the envelope at the later stage, the K-mechanism 
becomes to work efficiently, and leads to the stronger insta- 
bility. 

Although the g2-mode is also excited by the e- 
mechanism, its instability appears mainly for the lower mass 
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Figure 5. Properties of models and tlie dipole gi-mode at dif- 
ferent evolutionary stages for tiie 1.2Mq star, (a) work integral 
normalized with the total oscillation energy in a whole star, (b) 
radiative luminosity perturbation normalized with the value at 
the photoshere. (c) squared temperature perturbation normal- 
ized with the peak value in the deep interior, marked with the 
open circle, (d) effective temperature dependence of nuclear re- 
actions, £t, and ^He mass fraction {A color version of this figure 
is available in the online journal). 
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Figure 6. Work integrals for the gi and g2-modcs for different 
mass stars with Z = 1 X 10^*. 



stars. Since the damping just outside the nuclear burning 
core is stronger compared with the gi-mode, the g2-mode is 
not destabilized for the higher mass stars (figure [6|. For the 
lower mass stars, on the other hand, such damping is weaker, 
and the g2-mode is excited mainly by the e-mechanism. 

The g2-mode instability appears also in the late stage 
of ~ 1.6Mq stars with Z — I x 10~^. In this case, although 
the convective core size is not substantial, the contribution 
of the CNO-cycle becomes to exceed that of the pp-chain, 
and hence /i-gradient around the nuclear burning core be- 
comes steeper. In this situation, gravity waves are trapped 
in the steep /^-gradient zone, and the e-mechanism works 
efficiently. 
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Figure 7. Boudaries of instability regions for the dipole gi and 
g2-modes on the T^ff— g plane with different metallicity. The evo- 
lutionary tracks are for Z = Ix 10~*. The solid parts of the tracks 
correspond to stages at which the convective envelope is confined 
to the layer above r/R = 0.99, while the dashed parts to those 
at which the convective envelope breaks in below r/R = 0.99. 
The thick dashed line corresponds to the one in figure [2] {A color 
version of this figure is available in the online journal). 



4.2 Metallicity dependence 

Figure [7] shows boudaries of instability regions on the Teff- 
g plane for the gi- and g2-modes. They are cut off in the 
lower temperature side, where convective envelope becomes 
to extend below r/R — 0.99. In the Z = 1 x 10~^ case, the 
instability appears for stars with < 1.2Mq, and is induced 
mainly by the K-mechanism. As the metallicity decreases, 
the instability region, particularly for the gi-mode, extends 
toward the more massive stars. 

One of reasons for this is that metal-poorer stars have 
lower density contrast between the inner and outer regions 
than metal-rich stars. Figure |8] shows the ratio of the central 
density to the average density in a whole star, which repre- 
sents the density contrast, for Z — \ x 10""^ and 1 x 10"*. 
The ratio for metal-poorer stars is lower than for the metal- 
rich stars having the same mass, since metal-poorer stars are 
more compact and then have higher average density, while 
the central density is not substantially different among the 
metallicity. Due to this, the g mode amplitude is relatively 
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Figure 8. Ratio of central density to average density, pc/{p) for 
models with Z = 1 X 10~^ and 1 X 10~* at an evolutionary stage 
Xc ~ 0.55. The abscissa is stellar mass. 

larger in the deep interior, and the e-mechanism can work 
efficiently for metal-poorer stars. 

Another reason is that as the metallicity decreases the 
CNO-cycle becomes replaced by the pp-chain for more mas- 
sive stars. As discussed in Sec. 14.11 dominance of the pp- 
chain is favorable for the e-mechanism instability, while 
CNO-cycle burning keeps substantial size of a convective 
core during stellar evolution. 

Figure |§] shows regions corresponding to the instability 
with more than 50% contribution of the e-mechanism to 
the total excitation energy. Such instability appears with 
Z < 6x 10"*, and with log Tcff > 3.90 in the gi-mode case. In 
this temperature range, corresponding to the bluer outside 
of the classical instability strip, the K-mechanism of helium 
ionization does not work efficiently, and the e-mechanism is 
responsible for the instability. In the g2-mode case, on the 
other hand, such instability appears inside of the classical 
instability strip, while not in the bluer side except for the 
late stage of ~ I.QMq stars with Z = 1 x 10"*, at which 
the e-mechanism instability is induced because of the steep 
/i-gradient as described in Sec. 14.11 As shown in figure [51 
the g2-mode is strongly excited by the e-mechanism in the 
lower mass star case, while strong damping just outside of 
the nuclear burning core avoids the instability in the higher 
mass star case. 



5 DISCUSSION 

Due to the e-mechanism instability, the corresponding stars 
might exhibit pulsations with period ~ 1 hr. The g modes 
have relatively substantial amplitude at the surface as shown 
in figure \S\ although a nonlinear analysis is necessary to 
obt ain the absolut e value of the amplitude. 

iDziembowskl (|l982l . \l98t estimated the dipole gi- 
mode amplitude in the solar photosphere to be about 20 
cm s~^ by adopting t h e thre e mode coupling theory. In par- 
ticular, iDziembowskil l| 19831 ) considered the special cases of 
parametric resonance, in which a parent linearly unstable 
mode is coupled with two dauther linearly stable modes. 
The oscillation energy of the parent mode is converted into 
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Figure 9. Same as figure [T] but for the instability with more 
than 50% contribution of the e-mechanism to the total excitation 
energy. {A color version of this figure is available in the online 
journal) . 

that of the two daughter modes at a rate proportional to 
the product of the amplitudes of the three modes. The most 
effective for low-degree low-order parent g modes appears to 
be a coupling with a pair of similar high-degree daughter g 
modes to a state in which the amplitudes of all the three 
modes are steady. 

The value 20 cm s^^ is too small that the g mode 
oscillations are detectable for stars other than the Sun. 
On the other hand, SPB stars and 7 Doradus stars ex- 
hibit g mode osci ll ations with much higher amplitude. 
lAppourchaux et al.l l|201(]| ) suggested that such stars have 
only shallow convection zones, and therefore daughter g 
modes propagate much higher in the envelope and dissi- 
pate much more strongly, thus being limited themselves to 
much lower amplitudes and thereby being less able to ex- 
tract energy from their parent. Such a situation also might 
be valid for the g modes in the metal-poor low-mass main- 
sequence stars analyzed in this study, although the linear 
growth rate is much smaller than in the SPB and 7 Doradus 
cases. It is worth searching stellar pulsations induced by the 
e-mechanism since such pulsations have not been detected 
yet and would be a great tool to examine the first derivative 
of the nuclear energy generation rate. 

Metal-poor stars have attracted attention as clues to 
chemical evolution of the Universe. Thanks to advance in 
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the observational side , a lot of metal-p oor stars have been 
detected to date fe.g. ISuda et al.ll2008l ). Although most of 
metal-poor stars are very far from the Sun, faint metal- 
poor pulsating stars have become detected in addition to the 
Cepheids and RR Lyr stars in the other galaxies, e.g., 20 S 
Set stars at V ~ 23 in the Cari na dwarf spheroidal galaxy at 
~ 100 kpc (|Mateo et aLlliggsl ). Recently, higher number of 
faint variables are f ound in Large Mage llanic Cloud through 
the OGLE project (jPoleski et al.ll201(]| '). 

Most of metal-poor stars are thought to be consider- 
ably old or to have the equivalent age with the Universe. 
The e-mechanism instability reported in this paper, how- 
ever, appears in the early stage of the main-sequence. In 
some galaxies such as dwarf Irregulars, on the other hand, 
the chemical evolution might be delayed and metal-poor star 
formation regions remain still now. In such galaxies, young 
metal-poor stars may exist and are expected to be candi- 
dates for the e-mechanism pulsators. 

We also state a possibility that growth of the amplitude 
due to the vibrational instability might induce material mix- 
ing and have a significant influence on the later evolution of 
the star , as was expecte d conce rning the solar neutrino prob- 
lem bv iDilke fc Goughl (|l972l ). In particular, if the mixing 
occurred in the core, the surrounding cooler and hydrogen 
riched matter would be incorpolated into the core. This phe- 
nomenon might induce rejuvenation of the star and prolong 
the lifetime. To obtain a solution to this problem, we have 
to persue the nonlinear evolution. 



6 CONCLUSION 

We performed a linear nonadiabatic stability analysis of 
metal-poor low-mass main-sequence stars. We restricted our 
analysis to evolutionary models having negligibly thin con- 
vective envelopes, for which we can obtain a definite conclu- 
sion on the stability analysis without uncertainty in treat- 
ment of a convective envelope. We found that the dipole 
gi- and g2-modes certainly become unstable due to the e- 
mechanism for Z < 6 x 10"''. The mass range of this instabil- 
ity extends toward higher mass with decreasing metallicity. 
One of reasons for this is that for metal-poorer stars density 
contrast between the inner and outer regions is lower, and 
that the g mode amplitude is relatively larger in the nu- 
clear burning core. Another reason is that as the metallicity 
decreases the CNO-cycle burning becomes replaced by the 
pp-chain burning, which is favorable for the e-mechanism 
instability. As a result of the instability, stellar pulsation or 
material mixing is expected to occur. To confirm those, we 
require highly precise observations of faint metal-poor stars, 
and a nonlinear analysis of the oscillations. 
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